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We theoretically investigate the looping dynamics of a linear polymer immersed 
in a viscoelastic fluid. The dynamics of the chain is governed by a Rouse model 
with a fractional memory kernel recently proposed by Weber et al. (S. C. Weber, 
J. A. Theriot, and A. J. Spakowitz, Phys. Rev. E 82, 011913 (2010)). Using the 
Wilemski-Fixman (G. Wilemski and M. Fixman, J. Chem. Phys. 60, 866 (1974)) 
formalism we calculate the looping time for a chain in a viscoelastic fluid where 
the mean square displacement of the center of mass of the chain scales as t^^^. We 
observe that the looping time is faster for the chain in viscoelastic fluid than for a 
Rouse chain in Newtonian fluid up to a chain length and above this chain length the 
trend is reversed. Also no scaling of the looping time with the length of the chain 
seems to exist for the chain in viscoelastic fluid. 

I. INTRODUCTION 

Looping dynamics involving long chain molecules has been an active research topic in 
chemical physics [l-6|. Obviously loop formation is a primary ^tep in protein folding and 
thus finds lot of attention from chemists and biologists too 

ihe hydrodynamic interaction jo], [lo|, flexibility of the chain 

15| profoundly alters the looping time. But to our best knowledge the effect of viscoelastic 
fluid around a polymer chain on its looping dynamics has not been addressed yet. But 
it has been shown that the presence of viscoelastic fluid around the polymer changes the 
mode relaxation of polymers and thus it should also influence the looping dynamics. The 
importance of such study where a polymer is immersed in a viscoelastic fluid comes from 
recent experiments on biopolymers in viscoelastic environment such as cytoplasm [l6|. It 
has been shown very recently that the dynamics of the chromosomal loci in the viscoelastic 



7|, sL It has been shown that 
14j and the solvent quality 
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bacterial cells gets greatly affected and the diffusion becomes anomalous |l7l|. Diffusion of 
lipid granules inside viscoelastic cells have also found out to be anomalous [l8| . In this paper 
we address the problem of end to end loop formation for a chain immersed in a viscoelastic 
solvent where the chain dynamics is anomalous and is described by a Rouse model with 



fractional memory kernel 19|. For t 
polymer chain in a viscoelastic fluid 



lis we use the recently proposed model of a single 



19| . A Similar model have earlier been used in the 



20|. 



context of mode coupling theory of polymeric fluids 

The paper is arranged as follows. In section II the Wilemski-Fixman theory for the end to 
end loop formation is briefly discussed. The radial delta function sink used in the calculation 
is introduced in section III. Section IV deals with the model for a polymer in viscoelastic 
fluid. Section V presents the results and the paper ends with the conclusion in section VI. 



II. THEORY OF END TO END LOOP FORMATION: WILEMSKI-FIXMAN 

FORMALISM 

In this section we briefly discuss the Wilemski-Fixman formalism for the end to end loop 
formation in long chains. The dynamics of a single polymer chain with reactive end-groups 



is modeled by the following Smoluchowski equation 2l|, |22 |. 

^^^ffll^ = LP({R}, t) - kS{{K})P{{R}, t) (1) 

Here P{{R},t) is the distribution function for the chain that it has the conformation 

{R} = i?i,i?2, Rn at time t where Ri denotes the position of the ith monomer in the 

chain of n monomers. S{R) is called the sink function which actually models the reaction 
between the ends and thus usually is a function of end to end vector. L is a differential 
operator, defined as 



1=1 



d dU 



P({R},t) (2) 



Here Dq is the diffusion coefficient of the chain defines as the inverse of the friction 
coefficient per unit length and U is the potential energy of the chain. Wilemski and Fixman 
|l| then derived an approximate expression for the mean first passage time from Eq. ([T]). 
This mean first passage time is actually the loop closing time for the chain. Thus only valid 



3 



in the limit of infinite sink strength, /c — > oo 23|, |2J]. The expression for this loop closing 
time reads 



^Mt^-I) (3) 



\Cioo 

Where C{t) is the sink-sink correlation function defined as 

C{t) = J dnj rfRo^(R)G(R,t|Ro,0)^(Ro)P(Ro) (4) 

In the above expression G(R, t|Ro,0) is the the conditional probability that a chain with 
end-to-end distance Rq at time t = has the end-to-end distance R at time t; P{Ro) is the 
equilibrium distribution of the end to end distance since the chain is in equilibrium at time 

27j which depends only on the separation between the 



t = 0. S{R) is the sink function [25- 
chain ends. 

Now it is obvious that the knowledge of G(R, t|Ro,0) is prerequisite to calculate the 
sink-sink correlation function C{t) and hence the closing time r. In case of a flexible chain 
the Greens function and the end-to-end probability distribution functions are known and 
Gaussian. 

For a flexible chain the Greens function is given by 



3/2 



G'(R,t|Ro,0) 
Where 



2vr {R')^, 



(1 - <j>\t)) 



3/2 



X exp 



3(R-(/>(t)Ro) 



2(i?2) 



eq 



(5) 



(R(t).R(O)) 



eq 



(6) 



imeq 

is the normalized end-to-end vector correlation function for the chain. The above ensem- 
ble average is taken over the initial equilibrium distribution for end-to-end vector P(Ro). 

Similarly end-to-end equilibrium distribution for the flexible chain with Li = Nb^ (6 is 
the kuhn length and is the number of monomers) at time t = is given by 2l|, |28| 



27rLf 



3/2 



exp 



'2Lf 



(7) 
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With the above Gaussian functions the sink-sink correlation function can be written as 
a radial double integral. 

^-372 / ^'^R^S{R)dR / A'KRlS{R,)dR, x 

2txL^J (i_ (K^idY'^ Jo Jo 



exp 



3 + i?j 



2Li{l-ct>\t)) 



sinh [{?><p{t)RRo)/{Ll{l - 



(8) 



{m)RR^)/{Li{i-(^\t))) 

The above integral can be evaluated analytically for some specific choice of the sink 
unctions. With a radial delta function sink the above integral can be evaluated analytically 



29| 



III. THE RADIAL DELTA FUNCTION SINK 

For a radial delta function sink 



29| . S{R) = S{R — a), where a is the capture radius in 
the model the sink-sink correlation function defined in (Eq.([H])) can be evaluated exactly. 
Since in this case the integration over R and Rq can be carried out analytically, the looping 
time can be expressed in a closed form as follows. 



- ^ / exp [-(2xo0^(t))/(l - <p\t))] sinh [{2xo<P{t)) / {1 - <p\t))] ^ 



V 2a;o0(t)y^l-0'(t) 

with 

3a2 



2Ll 

Obviously if the end-to-end vector correlation function is known the closing time can 
be calculated by carrying out the integration over time. Throughout the paper it is assumed 
that the chain is in equilibrium at t = and only at t = 0"'' the viscoelasticity is turned on. 
Thus the equilibrium end-to-end distribution for the Rouse as well as for the Rouse chain 
with the fractional memory kernel is given by (Eq.Q). The last integration over time is not 
analytical so has to be carried out numerically. 

IV. MODEL OF POLYMER IN VISCOELASTIC FLUID 



It has been shown recently that a viscoelastic fluid background imparts a memory into 
the dynamics of the isolated chain immersed in the fluid. So the dynamics of a chain in 
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a viscoelastic fluid is modeled by a generalized Langevin equation 19|. In other words 
the dynamics of an isolated chain without self interaction and hydrodynamic interaction in 
a viscoelastic fluid is described by a Rouse model with fractional Langevin model kernel. 
Within this model, a linear chain described by a space curve R{n, t) has the following 
equation of motion 



n,t). 



(10) 



The random force /(n, t) is Gaussian |30| and satisfies the following fluctuation dissipation 
theorem jl^ 



(f (n, t)f (ni, ti)) = ^Kit - t,)6{n - n,)l. (11) 
where ^ is the friction coefficient and K{t — ti) is the memory kernel. 

This is exactly the model adopted by Weber et al. 19| to investigate the physics of subd- 
iffusive motion of a polymer in viscoelastic fluid. Although practically the same model has 
earlier been used by Min et al. 30|] in the context of fluctuation of the distance between a 



donor-acceptor pair within a single protein complex [31|, |32|. Also the Smoluchowski equa- 
tion (Eq.(IT])) no longer applies to the above chain, as the chain dynamics is non-Markovian 
in this case. 

In this case the normalized end to end vector time correlation function (Eq.(l6])) is given by 



0(t) 



(R(t).R(O)) 



eq 



eq 



16 (X,(t)X,(0)) 

p=odd 



eq 



(12) 



where the correlation function has been shown to have the following expression 19 1 



(X,(t)X,(0))^^ = 



N^T{3 - a 



(13) 



where kr, 



and Ea./3{x) is the generalized Mittag-Leffler function 33 1 



X 



i=0 



r(/3 + ai) 



(14) 



The generalized Mittag-Leffler function reduces to the regular Mittag-Leffler function 
for /3 = 1. In the limit a — 1, which is the case of Rouse chain in Newtonian fluid the 
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regular Mittag-Lefiler function becomes simple exponential and one gets back well known 
Rouse relaxation dynamics. But in a viscoelastic fluid for which < a < 1, (j){t) obviously 
decays non-exponentially and one would expect this to profoundly affect the end to end loop 
formation. In principle defined in Eq.( fT2i) should be put back into Eq.Q to carry out 
the integration over time to get the looping time. This would be the end to end looping 
time for a linear chain without any self or hydrodynamic interaction in a viscoelastic fluid. 
Unfortunately the regular Mittag-Leffler function (Eq. ffT^ with (3 = 1) can not be evaluated 
analytically for an arbitrary value of a. In general when a = 1/n where n = 2,3,4, .. the 
regular Mittag-Lefiler function can be expressed as a sum over incomplete gamma functions. 
We choose a = 1, which is the case for a Rouse chain in Newtonian fluid and a = 1/2 
which is a special case for a Rouse chain in a viscoelastic fluid the regular Mittag-Leffler 
function has analytically exact expressions. With a = 1/2 3^, the regular Mittag-Lefiler 

X 

function becomes £'1/2, i(x) = exp{x'^)Erfc{—x), where Erfc{x) = 1 — ^ / exp{—y^)dy is 



the complimentary error function. Using this exact expression for a = 1/2 one arrives at 
the analytically exact expression for 4>a=i/2i.'t)- 



K=i/2it)= Yl ^212^1/2,1 [-it/Tp,a=i/2y^^] = Yl exp (t/rp,„=i/2)^r/c( (t/rp,„=i/2) ^/^) 



p=odd p=odd 



where rp,«=i/^ ei^^^^- 



(15) 



Weber et al. 



19| showed that with a = 1/2 the mean square displacement of the center 



of mass of the polymer scales as t^/^. 

On the other hand for a Rouse chain in Newtonian fluid for which a = 1, the above 
expression becomes a sum over exponentials as all the odd modes relaxes exponentially. 



=iW = Y] ^-^exp(-t/rp,«=i) (16) 

p=odd 



where Tp,a=i = 



Notice 0Q,=i(t — )■ 0) = 1 as ^ = 1 and 4'a=i{t 00) = 0. Same is true for 

p=odd 

4>a=i/2 but it behaves as a stretched exponential at short t and inverse power law at long t. 
This can be seen by expanding exp{t / Tp^a=i/2) Er fc{{t / Tp^a=i/2Y^'^) at short t and doing an 
asymptotic expansion of the same at large t. 
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V. RESULTS 

The looping time with a radial delta function sink defined in (Eq.Q) is calculated for a 
chain in the viscoelastic fluid when the dynamics of the chain is described by a generalized 
Langevin equation (Eq. lfTOl) ). It is done as follows. First the end to end vector time corre- 
lation function is calculated for the chain using Eq. (fT5|) as only for this case with a = 1/2 
the correlation function can be evaluated exactly. Then 0q,=i/2(^) is put back in (Eq.Q) 
and the integration over time is carried out numerically to get the looping time. Plots of 
0a=i/2(^) (for the Rouse chain in Newtonian fluid) are shown in Fig. 1. Notice 

the same set of parameters are used along with the same value of the friction coefficient (^). 
As one can see that the initial decay of (pit) is faster in viscoelastic fiuid than in Newtonian 
but beyond a time scale the trend is reversed and 4'a=i/2{'t) approaches zero extremely slow. 
Also for a shorter chain the initial decay of 0^=1/2 (^) is even faster. Faster the decays 
faster the integrand in Eq.([9]). Thus for a short chain immersed in the viscoelastic fiuid 
the contribution from the initial decay is so small that it did not get compensated from the 
longtime contribution and results a faster looping. While for a long chain the longtime dy- 
namics contributes appreciably so that the area under the integrand is big enough to make 
loop formation slower. This can be seen from the plot of /n(r) vs ln{N) in Fig. 2. Another 
important observation is the absence of any well defined scaling of the looping time with 
the chain length for the chain in viscoelastic fiuid (a = 1/2). This is also evident from the 
/n(r) vs ln{N) plot in Fig. 2. But the Rouse chain in the Newtonian fiuid [a = 1) shows a 



T ~ A^^ scaling as expected [10 



15|. 



VI. CONCLUSIONS 



In this paper we calculate the looping time for a chain immersed in a viscoelastic fiuid 
where the mean square displacement of the center of mass of the chain scales as t^^^. We 
found that up to a chain length viscoelastic fiuid actually enhances the looping rate as 
compared to a Newtonian fiuid with the same friction. But beyond a chain length the trend 
is reversed. This observation is due to the faster short time decay of the normal modes of 
the chain in the viscoelastic fiuid as compared to in a Newtonian fiuid. This is eve more 
prominent for short chains. Also no scaling of the looping time with the chain length seems 
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FIG. 1: against time (t). The values of parameters used are: = 100, b = 0.1,^ = 1, ksT = 1. 

to exist in viscoelastic fluid. Here we would like to point out that in a real viscoelastic fluid 
a chain would definitely feel a larger frictional drag. Thus in a real viscoelastic fluid the 
chain would relax even slower and one would expect slower looping time as compared to 
normal Rouse chain even with short chains. 

In future we would like to investigate the effect of self interaction along with the vis- 
coelastic background on the looping dynamics. It is obvious that in a crowded environment 
self interactions will also have profound influence on the looping dynamics of the chains. 
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FIG. 2: /n(r) against ln{N). The values of parameters used are: a = 0.5,6 = 0.1,^ = 1, fesT = 1. 
project of K. L. Sebastian. The author thanks K. L. Sebastian for encouragements. 
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